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ABSTRACT 

At intervals as short as ten thousand years, each white dwarf (WD) passes 
within a solar radius of a planetoid, i.e., a comet, asteroid, or planet. Gravi¬ 
tational tidal forces tear the planetoid apart; its metal-rich debris falls onto the 
WD, enriching the atmosphere. A third of WDs exhibit atmospheric “pollution”. 
For roughly every hundred planetoid disruptions, a planetoid collides with a WD. 
We simulate a small number of collisions, in which “death-by-dynamics” refers to 
the fate of the planetoid. We also compute the energies and likely durations of a 
broad sample of collision events, and identify detection strategies at optical and 
X-ray wavelengths. Collisions with the most massive planetoids can be detected 
in external galaxies. Some may trigger nuclear burning. If one in 10^ - 10® of 
WD-planetoid collisions creates the conditions needed for a Type la supernova 
(SN la), “death-by-dynamics” would also refer to the fate of the WD, and could 
provide a novel channel for the production of SN la. We consider the circum¬ 
stances under which the rate of SNe la can be increased by interactions with 
planetoids. 


1. Planetoids, Collisions with White Dwarfs, and the Generation of Type la 

Snpernovae 

1.1. White Dwarfs and Tidal Disruptions 

Planets around white dwarfs (WDs) have yet to be discovered, but several lines of 
reasoning suggest that WDs commonly host planetary systems (Di Stefano 2011; Parriott 
& Alcock 1998). The strongest evidence is that more than 1/3 of WDs exhibit substantial 
atmospheric “pollution”. The mechanism responsible appears to be the tidal disruption 
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(TD) of planetoids (Jura 2008; Zuckerman et al. 2010; Farihi et al. 2013). 

Because the metal settling time can be ~ 10^ years (Koester 2009), contributions from 
planetoids can be assessed only in aggregate. The average mass influx is estimated to 
be iif = (10^^ — 10^^)g/yr (Zuckerman et al. 2010; Barber et al. 2012). Some WDs may 
therefore accrete up to an Earth mass of metal-rich material during a Hubble time. 

The mass is deposited by sequences of incident planetoids: comets, asteroids, and 
planets. To relate M to the influx of planetoids, we use a power-law mass distribution, 
oc M~^, with /3 = 1.75. Here, dN(M) is the number of objects with mass between 
M and M dM. Thus, the measured rate of pollution can be provided by the disruption 
of one 10^^ g TD per year, one 10^^ g TD per century, or one 10^^ g TD per Myr. The mass 
distribution predicts a mix of masses and time intervals between events. 


1.2. Collisions 

If planetoids come close enough to WDs for TDs, then some must collide with WDs. 
Because the cross section for close approaches is proportional to the distance of closest 
approach, the collision rate is ~ 1% the rate of TDs, scaling roughly as the ratio Rwd/Rq, 
where Rwd is the radius of the WD, and Rq is very close to the tidal radius at which 
disruption occurs (e.g., Pineault & Landry 1994; Schlichting et al. 2013; Puentes et al. 
2010). The word “death” in our title can therefore refer to the certain destruction of a 
planetoid when it collides with a WD. In §2 we show the results of a small number of 
simulated collisions and then compute the detectable effects of a larger number across the 
mass range spanning from comets to planets. 
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1.3. Type la Supernovae 

In extreme cases, a sequence of TDs and collisions may spark nuclear burning, 
potentially providing a novel pathway to SNe la. If interactions with planetoids can help 
propel WDs to SN la explosions, the word “death” in our title may also refer to the end 
of the WD. The addition of a new evoluionary pathway may be a positive step toward 
developing a better understanding of how SNe la are generated. Observations tell us that 
of every hfty WDs, one ends as a Type la supernova (SN la). Yet calculations based on the 
standard models have difficulty producing rates this high (Nelemans et ah 2013). 


2. Collisions between WDs and Planetoids 
2.1. Energy of WD-Planetoid Collisions 

An incoming planetoid is subject to strong tidal forces once it crosses its tidal radius 
rt = RwDiPwD/PpiY^^, where pwo, Ppi are the average densities of the WD and planetoid, 
respectively. At this moment, the planetoid disintegrates (likely into chunks) and has 
a spread in binding energies AE = G that translates to a spread in times 
of impact At ~ (l?pi/'rt)tff, where is the amount of time elapsed between when the 
planetoid crosses rt (~ 0.5 i?© for Mwd = Af©) and when it impacts the WD’s surface, 
iff = TT rf/SGM wd 10 min. For a 4 km planetoid with ppi = 10 g cm^ this spread would 
only be 7 ms, corresponding to an impactor whose length has increased by a factor of ~ 5. 

This modest increase in length means that the impactor strikes the surface of the WD 
in a small area, comparable to its own size, rather than being spread over the WD’s surface 
a la the Shoemaker-Levy comet which impacted Jupiter on its second periapse (Asphaug & 
Benz 1996). We simulate a two examples of planetoid impacts using the hydrodynamical 
code FLASH (Fryxell et al. 2000), the results of which are presented in Figure 1. We hnd 
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that the planetoid penetrates to a depth where the density of the WD’s atmosphere is a 
factor of ~ 10 greater than ppi, with the snrface and the impactor being heated to a few 
million degrees. After impact, the ejecta qnickly expands and cools, initially radiating 
in the X-rays for ~ tens of milliseconds for the 4 km impactor, longer for more massive 
impactors, bnt releasing the majority of its energy in the optical over longer timescales. 


2.2. Event Detection 

State-of-the art optical snrveys - e.g., Pan-STARRS (Kaiser et ah 2002), Palomar 
Transient Factory (PTF; Ran et ah 2009), SkyMapper (Keller et ah 2007), GAIA (Perryman 
et ah 2001), Palomar-Qnest (Djorgovski et ah 2008), Catalina Real-Time Transient Survey 
(Drake et ah 2009) - can reach a held depth up to roughly ~21 mag per observation 
(in one-to-several minutes exposure time; Ran et ah 2009, and references therein). We 
adopt the combined characteristics of Pan-STARRS and PTF as our benchmark. Current 
surveys yield sky coverage up to several thousands of square degrees each night. To 
classify a collision event as robustly detectable by present surveys, we adopt a strong 
threshold of 20 mag, in order that non-detections in data bracketing the transient event 
can be discriminating. To estimate an event’s detectability, we adopt a time of 100 s for a 
ground-based exposure and neglect interstellar extinction. We assume 100% efficiency in 
the conversion of gravitational potential energy into optical radiation (for simplicity, we 
neglect any bolometric correction). The emission timescale is calculated following Arnett’s 
law (Arnett 1980) for an expanding gas. Recombination is the dominant source of emission, 
and continues until the gas has fully recombined, a timeframe which scales linearly with the 
size of the impactor, with a rule of thumb ~ 1 s/km. The effective optical luminosity Lopt 
is the luminosity an observer would infer from a 100s exposure around the peak. Note that 
when the emission timescale becomes very short, the effective luminosity will be lower than 
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Fig. 1.— Results from two hydrodynamical simulations of the impact of a tidally stretched 
bolide with a WD atmosphere. The left three panels show 2D slices of log p as a function of 
time t through the plane of symmetry for an impactor initially 400 km in radius, where the 
impactor travels right to left at a 45° angle to the WD’s surface. The middle three panels 
show the same for a 4 km impactor. The right panel shows a 3D snapshot from the 4 km 
impactor simulation, with the contour corresponding to the region that contains at least 10% 
of the impactor’s composition by mass, color-coded by temperature. 
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the short-lived peak luminosity. 

While optical surveys cover regions of the sky at a specihed cadence, X-ray monitors 
(e.g., Swift BAT, Barthelmy et al. 2005; RXTE ASM, Levine et ah 1996; MAXI GSC, 
Mihara et al. 2011; Integral IBIS, Ubertini et al. 2003), often cover nearly the full 
sky continuously, producing ~daily maps of the brightest X-ray sources. One of the 
most-sensitive of these monitors, the MAXI GSC, can detect source fluences of ~ 50 
Crab seconds per orbit (i.e., ~ 10“® erg cm“^, corresponding a persistent X-ray flux < 50 
mCrab ~ 10“® erg s“^ cm“^). MAXI observes very nearly the full sky each ~ 90 minute 
orbit. We adopt this fluence as a detection threshold. (Because the prompt emission is 
short compared to the exposure timescale, fluence and not flux is the critical quantity.) 

In calculating the X-ray emission, we arbitrarily assume that a fraction fx = 10% of the 
impact energy E is released in the form of prompt (nonthermal) X-rays. 

By contrast, Chandra / XMM-Newton imaging is sensitive to much fainter X-ray 
sources. The tradeoff, of course, is that these instruments have very narrow £elds-of-view, 
and a correspondingly low duty cycle for observing any particular region of sky. We 
neglect photon pileup effects, and adopt a detection threshold of 25 X-ray photons (at a 
characteristic energy of 1 keV, for an instrument with effective area ~ 300 cm^). 

In Table 1, we compute rough estimates for the rate, brightness, and detectability of 
impact events across a spectrum of masses. The impactor masses we consider range from 
10^^ g, based on a lower optical threshold of Lopt = 10^® erg/s, to Jupiter’s mass, 10^° g. 
In each case, we give the maximum distance for which the event could be detected in the 
optical and X-ray detectors discussed above. 

Asteroids of mass < 10^°g produce the least energetic events in Table 1. These 
are most readily detected in the optical, just within our Galaxy (or even just the Solar 
neighborhood). Although event rates are high, so that hundreds to thousands of events are 


Table 1. Event Rate and Detectability 
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detectably bright, these transient impacts are so brief that even high-cadence ground-based 
surveys designed to hunt transients are only likely to hnd 0.1-1 per decade (being most 
sensitive to impactors with M lO^^g). 

Meanwhile, energetic events involving sub-Earth-mass planets (M ~ 10^'^g) occur 
with sufficient frequency that they are potentially detectable with X-ray imaging of galaxy 
clusters. In particular, several Ms aggregate imaging is available for crowded helds of 
nearby rich clusters {d < 100 Mpc, e.g., Virgo, Norma, and Coma). Our model predicts 
<0.1 events in the aggregate data set. However, we speculate that the vigorous dynamical 
activity in clusters may enhance these rates. 

Intermediate impactor masses just below the scale of Ceres (M < 10^"^ g), are 
sufficiently faint that they are all but undetectable outside the Local Group. At the same 
time, the collision rates are so low that pointed ground or space-based instruments have 
negligible chance at serendipitous discovery. In this regime, the X-ray monitors shine. 
Nearly 1 event per year is detectably bright by an X-ray monitor; however, factoring in the 
sky coverage, only 0.01 — 0.1/yr detections are expected. 

Next decade, LSST, with its tenfold improvement in sensitivity is well-poised to 
make dozens of detections of events by both the lowest-mass and highest-mass mass 
impactors. In particular, we predict that tens of Galactic events involving impactors with 
M = (10^^ — 10^®) g will be detected in the hrst year. The dominant sensitivity is for 
masses M ~ 10^® g, which roughly corresponds to the lowest-mass impactor for which the 
detectable volume is matched to the scale of the Galaxy. We predict tens of planet impact 
detections each year: we expect an annual rate of roughly one Mercury-scale detection, and 
nearly ten Jupiter-mass event detections. In the latter instance, the signal is dominated by 
high redshift events (;2 ~ 1 — 3). 


In short, we hnd that with current instruments, WD-planetesimal impacts are at 
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the cusp of detectability. It is unlikely that any event has heretofore been identihed or 
recorded. At the same time, we note that with a generational improvement in sensitivity 
for X-ray coded-mask instruments, or in the forthcoming LSST at optical wavelengths, we 
will be well-poised to discover such impact events. Notably, a modest improvement in X-ray 
monitor sensitivity would critically extend the detectability of Jupiter-WD collisions out to 
distances of > 100 Mpc, the proximity of nearby galaxy clusters. Detection rates for these 
most energetic events would then become appreciable. Even as we wait for a statistical 
sample of these events to be established using the rich capabilities of LSST, in this decade 
monitoring programs by PTF, Pan-STARRS, and other campaigns such as OGLE will be 
increasingly capable of discovering a nearby planetoid-WD impact. 


3. The Progenitors of Type la Supernovae 
3.1. Possibilities and Caveats 

WD-planetoid collisions are frequent enough that if only one in ~ 10^ — 10® produces 
an SN la, death-by-dynamics models will contribute signihcantly to the total rate. To 
trigger an SN la, a collision must dynamically ignite matter on a WD’s surface. This 
requires compressing and heating material above a critical temperature where the burning 
time becomes comparable to the dynamical time. For pure hydrogen, the critical mass to 
auto-ignite a thin atmosphere is small, and such WDs have very little in their atmospheres 
(no more than 10“^ Mq). For pure helium, the maximum envelope mass can be tenths of 
a solar mass for low-mass WDs (~ O.SM©), and decreases with increasing accretor mass 
(Bildsten et al. 2007). 

If, however, a part of the envelope is compressed and/or heated above its steady-state 
condition, the fluid can spawn a self-supporting nuclear detonation that propagates 


11 


throughout the envelope. Initiating such a detonation involves heating a mass far smaller 
than the total envelope mass. For pure helium, this mass is on the order of Earth’s mass. 
Hence, we would expect that, at minimum, the total impacting mass required to initiate a 
detonation in the helium layer should be equal to this value. Furthermore, the conditions 
for detonation are usually realized at the interface between the outer envelope and the 
inner core. In cases where the outer envelope mass is large, an impactor would need to 
penetrate through a mass greater than its own to reach the depth of this interface. Secondly, 
planetary/planetesimal impactors are low-density, an Earth-mass body has roughly the 
same size as the WD, and its kinetic energy is spread over a large shell (meanwhile, the 
critical length scales are derived presuming all the mass lies within a small, spherical 
volume). 

We envision two circumstances which could change this picture. Firstly, the critical 
length scales have only been calculated for a few composition combinations, namely 
carbon-oxygen (Seitenzahl et ah 2009), and pure helium (Holcomb et al. 2013). However, 
the inclusion of even a small admixture of heavy nuclei into a helium layer will greatly 
reduce the burning timescale, since a captures onto these nuclei are far more rapid than the 
relatively slow 3-a reaction (Shen & Moore 2014). Consequently, while a single impact may 
not detonate a pure helium envelope, repeated bombardment may contaminate the envelope 
to such a degree that a captures and other nuclear pathways become viable. Secondly, 
because the impact is fundamentally a dynamical process, the geometry of the impact may 
lead to a shock-focused detonation in a layer stable to a single hot spot generated by the 
collision - see e.g., Fink et al. (2010). 
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3.2. Death by Dynamics in Tandem with SD and DD models 

The two leading models for SN la progenitors invoke interacting binaries with epochs 
of mass transfer and common envelope evolution. These epochs last long enough to permit 
sequences of planetoid interactions. In single degenerate (SD) models, a non-degenerate 
companion donates matter to the WD. In double degenerate (DD) models two WDs merge, 
leading to a detonation. 

Planetoid interactions that enrich the WD may alter SN la rates in three ways. Firstly, 
if enrichment extends (narrows) the range of infall rates compatible with steady nuclear 
burning (Piersanti et ah 2014), then both DD and SD rates would be affected. Secondly, 
enhanced metallicity increases winds, possibly stabilizing mass transfer for systems verging 
on dynamical instability. Finally, enrichment may make it easier or harder to spark 
explosions during the accretion/merger processes associated with explosions in SD/DD 
models. Furthermore, a subset of planetoid interactions with the donor, disk, or common 
envelope could alter the binary’s further evolution. 


3.3. Predictions of Death-by-dynamics Models 

Death-by-dynamics models predict diversity among SNe la. The WD mass does not 
have a hxed value at explosion. Differences in planetoid compositions produce chemical 
diversity. These models are novel in positing that a WD can become an SN la without 
a close companion. There may, however, be other stars near the WD; distant bound 
companions, passing stars which have triggered planetoid interactions, or even close 
interacting companions. Thus, the properties of any stars in the held close to the SN 
may exhibit a variety of luminosities and spectral types. The circumstellar region may be 
enriched with planetary debris that interacts with the supernova. 
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Early-time explosions are predicted by death-by-dynamics models, because stellar 
evolution destabilizes planetoid orbits. As well, planetoids continue to be placed on 
eccentric orbits over a Hubble time. The combination of strong early-time signals and 
continuation through late times is consistent with data on SNe la (e.g., Sullivan et al. 
(2006)). Calculations are needed to check whether the distribution of death-by-dynamics 
delay times is consistent with observations. 

The discovery of WD-planetoid collision events will identify environments that produce 
SNe la through death-by-dynamics channels. Conversely, regions producing SNe la must be 
rich in WD-planetoid interactions. 


4. Conclusions 

WD-planetoid collisions produce bright, detectable events. Observing these events will 
measure peak luminosities and time evolution in broad energy bands, and will quantify 
environmental effects. This will provide clues to the mass distributions of both planetoids 
and WDs, and eventually chemical information. Such studies are possible because large 
numbers of planetoids orbit WDs and large numbers of WDs inhabit galaxies. Even unusual 
interactions and sequences of interactions are common in galaxies and galaxy clusters. 

A planetoid interaction or sequence of them may spark SNe la. The exploding WDs 
may be isolated or may be in interacting binaries. While it will be challenging to test 
death-by-dynamics SNe la progenitor models, we note that this has also been true for SD 
and DD models. Fortunately, if death-by-dynamics models do produce SNe la, they make 
a host of predictions, some of which can be tested through the study of bright events 
caused by WD-planetoid collisions. In summary, death-by-dynamics models make testable 
predictions and suggest intriguing directions for future investigations. 
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